Dynamical Transitions of a Driven Ising Interface 
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We study the structure of an interface in a three dimensional Ising system created by an external 
non-uniform field H{r,t). H changes sign over a two dimensional plane of arbitrary orientation. 
When the field is pulled with velocity v^, (i.e. H{r,t) = i/(r — Vei)), the interface undergoes a 
several dynamical transitions. For low velocities it is pinned by the field profile and moves along with 
it, the distribution of local slopes undergoing a series of commensurate-incommensurate transitions. 
For large Ve the interface de-pinns and grows with KPZ exponents. 

PACS numbers: 68.35.Ja, 05.10.Gg, 64.60.Ht, 68.35.Rh 



Introduction: The study of the structure and proper- 
ties of surfaces and interfaces in condensed matter is of 
considerable technological importance^, Many phase 
transitions and chemical reactions originate at a surface 
and their progress in time is often synonymous with the 
motion of resulting parent- product interfaces Q. Often, 
interfaces undergo structural transitions, like reconstruc- 
tions and roughening, the study of which have recieved 
considerable attention in the past because of its impor- 
tance to a variety ofprocesses like catalysis and surface 
phase transitionsP, 3]. Surface structure can also be in- 
fluenced by its dynamics. Consider a simple interface be- 
tween the ground states of an Ising model[5] with nearest 
neighbor ferromagnetic interactions defined on the d di- 
mensional cubic lattice. The equilibrium, static, interface 
defined by the height h{r) fluctuates so that the mean 
squared deviation of h from the average grows @ as 
and t'^ where L is the typical system size, t the time with 
the exponents a = (2 — d)/2 and /3 = 2 — d/4. Above 
d = 2 the interface remains flat {a ~ f3 — 0). On the 
other hand, a moving interface in this system when driven 
by an uniform field coarsens such that a ^ 2/{d + 3) and 
/3^1/id + 2). 

Motivation: In a few earlier publications 0, @, we 
have shown that it is, nevertheless, possible to stabilize 
a macroscopically flat (a = (3 = 0) but moving interface 
in this system by using a non-uniform field (see Fig[T]) 
of the form H{r,t) = i?otanh[(r — v^t) ■ n/x] where the 
unit vector n fixes the orientation of the interface and x 
the (intrinsic) width. In Ref.[l, 0] the properties of such 
an interface in the two dimensional Ising model was ana- 
lyzed in detail. The main results of this study were as fol- 
lows. Firstly, there exists two distinct dynamical phases: 
for small Ve, the interface, defined as the locus of all 
the points where the magnetization S{r, t) changes sign, 
moves together with the field with velocity Vf = Ve, and 
its average position is fixed at all times by the plane {line 
in two dimensions) ZQ{x,y) over which the field changes 
sign and the interface remains macroscopically flat. For 
large Ve > foo, on the other hand, the interface detaches 
from zq and grows with velocity Wqo essentially driven by 



FIG. 1: (left) The geometry of our simulation cell. Ising spins 
are defined over a 14 x 14 x 30 lattice elongated in the z di- 
rection. An external field is imposed which changes sign over 
a plane A, B, C making angles 6 and (j) with the x and y axes 
respectively. The plane A, B, C moves to the right with veloc- 
ity We. (right) A plot of the velocity of the Ising interface vj 
as a function of Ve for two different orientations - (0,2/7) (top) 
and (2/7,2/7) . Note that for small Ve, the interface follows 
the field {vf — Ve) but gets detached at higher velocities. 



an uniform field at the same time coarsening with KPZ 
exponents The value of Voo is orientation depen- 

dent because the motion of interfaces at low temperatures 
takes place mainly by flipping the relatively high energy 
spins occupying corner sites. The interfacial orientation 
with the largest density of such unstable spins has the 
highest Woo- 

Secondly, the low velocity pinned phase was shown to 
be rather interesting because it shows an infinite number 
of structural phases and phase transitions very similar 
to commensurate-incommensurate transitions in adsor- 
bates on solid surfaces 11, 12|. In the pinned phase, 
the spatial profile of the field forces the magnetization 
to follow H as closely as possible. Specifically, the av- 
erage orientation of the interface is fixed by the orienta- 
tion of the line zq over which H changes sign. Since the 
underlying lattice admits only a discrete set of rational 
orientations, at equilibrium, the interface settles down to 
the closest rational approximant possible given the finite 
size of the system. The pinning energy depends on the 
orientation of the interface relative to the lattice, simple 
rational fractions being more strongly pinned. Increas- 
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ing Ve has the effect of increasing dynamical noise, and 
orientations with stronger pinning and shorter relaxation 
times become preffered. As is increased, therefore, an 
interface with an arbitrary orientation distorts locally so 
that over most of its length, the orientation conforms 
to a low order rational fraction. A non-zero density of 
rfisco77i77ie7isMra<zons0, 0] maintains the average orien- 
tation of the interface equal to that of zq- A plot of 
the most probable orientation of the interface against 
the average orientation for any Ve < v^o therefore shows 
steps correspond ing to an incomplete Devil's staircase 
structured, 

The purpose of the present Brief Report, is to deter- 
mine whether these conclusions are valid for the d = 2 
interface in the c? = 3 Ising system and are not an arti- 
fact of the reduced dimensionality. Most experimentally 



relevant interfaces being two dimensional, we believe this 
to be an essential point which needs to be checked by di- 
rect calculation. We show that the general scenario per- 
sists even in the higher dimensional case studied here, 
although there are some special features related to the 
larger degree of freedom of the two dimensional inter- 
face compared to the one dimensional case studied in 
Refs.[l,l5- Unlike the former study, where one may make 
use of a formal mapping of the Ising interface dynamics 
problem to an Assymetric Exclusion Process [ISf to de- 
velop efficient numerical codes and obtain useful analyt- 
ical results, no such mappng is possible for the higher 
dimensional case and we use standard Monte Carlo sim- 
ulations using Metropolis update for this problem. The 
exact nature of the update rule is not expected to inffu- 
ence the qualitative nature of our results. 




FIG. 2: (left (a) and middle (b)) Plots of the averaged interface defined by the locus of zeroes of the magnetization for two 
different values of Ve = .1 (left) and .0125. The orientation of the interface is (2/7,2/7). Note that for small smaller Ve the 
interface is more sharply defined, (right (c) ) A portion of the dynamical phase diagram for a one dimensional interface taken 
from Ref.^^] showing locked (white regions), fluctuating (black regions) and de-pinned (hatched region) phases. Note that an 
interface with orientation close to (and to the right of) 2/7 would pass through regions where the most probable orientation 
shifts to (1/4) as the We is increased. 



Mean field results: When spatial flucutations of the 
interface are neglected, the velocity of the interface w/, 
defined as the locus of all the points where the magne- 
tization ^(r, t) changes sign, is given by the solution of 
'Vf — fooSign(ti/ ~ Ve) at late times Q. There are two 
regimes: (a) for small Ve < Voc the interface is pinned 
by the field and its position coincides with the plane over 
which the field changes sign so that Vf = and (b) for 
Ve > Voo, the interface depins and grows with its orien- 
tation dependent intrinsic velocity Voo ■ These mean field 
results are independent of dimensionality and should hold 
for Ising systems in any dimensions. We show, indeed, 
that this is true and the depinning transition, predicted 



by the mean field calculation persists in 3 — d (see Fig[T]) . 

Simulation details: Owe system consists of Nx{— 14) x 
Ny{= 14) X N^{= 30) Ising spins arranged in a simple 
cubic lattice [14| interacting ferromagnetically with their 
nearest neighbors. The boundary conditions for this sys- 
tem has to be specified with some care since we need to 
describe, in general, oblique interfaces with arbitrary ori- 
entations. We are interested in the low temperature, high 
field limit, however, which greatly simplifies this task. In 
this limit, only a few layers of spins immediately adjacent 
to the interface have any nontrivial dynamics, spins far 
away from the interface being frozen in the direction of 
the local external field H. We use shifted periodic bound- 
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ary conditions ensuring that the periodic image of each 
spin has the same value of H (see Fig 1.). In the z direc- 
tion we have open boundary conditions the exact nature 
of which is irrelevant as the spins at the two extremeties 
are, in any case, frozen to the values fixed by the sign of 
H. 

A typical run consists of equilibrating the spins in a 
given updating H according to the value of Ve mea- 
sured in units of lattice parameter per Monte Carlo step 
(MCS) and computing the non-uniform magnetization. 
The value of the magnetization is stored for each time 
step till the interface traverses the length of the sample. 
This process is typically repeated 10"^ times and the av- 
eraged time dependent magnetization profile is obtained. 
To obtain the location of the interface the averaged mag- 
netization S{r, t) is fitted to a general hyperbolic tangent 
profile. The locus of zeroes is the interface. The proba- 
bility distribution of slopes of this interface is then com- 
puted using finite differences and is further averaged over 
10^ independent runs. Therefore, although our system 
size of 5880 spins is modest, the scale of our computa- 
tions is not; a run with a single parameter value needs to 
be repeated many times to gather sufficient statistics. 

Results and Discussion:We begin our discussion of re- 
sults by considering first the velocity of the interface v / 
as a function of Ve- We obtain this from the slope of a 
plot of the position of the center of mass of the inter- 
face as a function of time. Our results for two different 
orientations of n are shown in Fig. 1 (right panel) . The 
depinning transition and saturation of to a orientation 
dependent value of Voo is clearly visible. The saturation 
value Voo depends on the density of steps which, in turn, 
is a function of the interfacial orientation. Compared 
to the two dimensional case of Refs.jl, [sl, however, this 
transition appears to be less sharp. This is, in fact, a 
finite size effect d, 14 1; the sharp transition is recovered 



only in the large time limit which is difficult to achieve 
for our modest system size. Note that the transition is 
broader for the faster moving interface for which the ob- 
servation time is shorter. 

We concentrate next on the interfacial structure of the 
pinned interface. The average orientation of the interface 
is given by the unit vector fi which may be parametrized 
by the angles 9 and (f> as shown in Fig. 1 (left panel). We 
find it more convenient to define the slopes me = tan(6') 
and = tan(<^). The structure of the interface aver- 
aged over many runs is shown in Figl2] (a) and (b) for 
two different velocities and for riig — — 2/7. It is 
clear from these pictures that the velocity strongly in- 
fluences the structure, smoothening it on the average. 
In order to understand better this smoothening process, 
we compute the probability distribution of the local and 
instantaneous slope using a resolution window of 5 x 5 
spins. Within this window, the local slope is determined 
at every intance of time and the data is used to build up 
the joint probability distribution P{mg, m^). 



The joint probability distribution P{me, m^) as well as 
the projected distributions P{mg) and P{mcj,) are plotted 
in Figl3]for three different Ue- In order to understand the 
sequence of the transitions, we refer to a portion of the 
dynamical phase diagram in the Ve — me plane Fig|5]Jc). 
In this figure, the white regions represents the lock-in 
phases where the most probable orientation locks in to 
the closest rational slope which may be different from the 
average slope of the interface. The black regions repre- 
sent incomensurate states where the the most probable 
slope is nearly equal to the average slope and the inter- 
face dynamically fiuctuates between nearby orientations. 
Finally, the gray region for large corresponds to the 
depinned interface. An interface with a slope me close to 
but smaller than 2/7 would first lock-in at 2/7 for small 
Ve and then to 1/4 and finally depinn as Ve increases. 
In-between these transitions the interface would be in- 
commensurate . 

This sequence of transitions is apparent in the plots of 
P{me,m^) shown in FigO However, in contrast to the 
two dimensional situation study Q, the interface in the 
three dmensional Ising system has more degrees of free- 
dom. This is evidenced by the presence of multiple local 
orientations which show up as secondary peaks for small 
Ve . Careful observation of the structure of the interface 
shows that the local slope of the interface continually 
fiuctuates in space and time over a small region of phase 
space around the average orientation. 

What determines the strength of the pinning in any 
particular orientation? The residence time in any par- 
ticular orientation depends both on the strength of the 
pinning as well as on the nearness in configuration space 
from the average orientation. In the study of the one di- 
mensional interface, the pinning energies could be com- 
puted analytically by mapping the one dimensional in- 
terface to a sequence of particles in a one dimensional 
lattice[l3, Ell- This procedure is not possible in the 
present case. However, the conclusions of this study 
seems to be valid for each of the two slopes me and 
7710 taken separately. A simple ansatz which proceeds 
to break up the response of the 2 — d interface with slope 
(me,m^) as that of two 1 — d interfaces with slopes me 
and m^ succeeds fairly well. The lowest energy orien- 
tations correspond to 777^ or m-^ being of the form 1 /n 
where n is an integer. So the interface with an initial 
slope near to 2/7,2/7 (FiglSjc)) for the smallest veloc- 
ity, breaks up into regions with slopes close to 2/7 with 
the largest probability corresponding to a slope close to 
1/4, 1/4 (Fig. [3)(b)) which is the similar to that for the 
1 — d interfacefSj. However, as mentioned before and as 
is clear from Fig. ^h). there are considerable weights 
for nearby slopes where the interface spends a consider- 
able amount of time. These appear as secondary peaks 
in P[me,m^). If the velocity is increased further, the 
interface becomes incommensurate and the probability 
distribution nearly Gaussian with a sharp peak only for 
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FIG. 3: Projected and two dimensional distributions of the 
local slope of the interface with an average slope close to 
mg = 2/7 and = 2/7 for various values of Ve- Left panel 
shows the probability distribution of the slope P{me), the 
middle panel P{m^) and the right panel the full two dimen- 
sional distribution P(mg, m^) for Ve = .2 (a), .1 (b), and .0125 
(c). Note that in (a), the most probable slope is close to the 
average slope of the interface. As the velocity is lowered, the 
interface tends to lock in to a rational slope close to 1/4, 1/4, 
for which the probability shows the strongest peak (b). As the 
Ve is lowered further, the lock-in slope is a higher order frac- 
tion 2/7, 2/7 which is closer to the average orientation of the 
interface. Large fluctuations away from the lock-in slope are, 
however, present throughout contributing secondary peaks in 
the probability distribution function. 



the average slope (Fig. [S^a)). Further we have verified 
with actual computation that the 2-d probability distri- 
bution is not given by a product of the 1-d probabilities 
making fluctuations in mg and highly correlated. 

A detailed analytical calculation of the full probability 
distribution would be useful but is beyond the scope of 



this work. Future work in this direction would hopefully 
explain these issues analytically and obtain the full three 
dimensional phase diagram in the mg,m(j) and plane. 
We hope our studies will be useful in understanding mov- 
ing interfaces in more realistic systems 16| and t he g rowth 
of colloidal crystals using patterned substrates [17 1. 

The authors thank A. Chaudhuri, A. Milchev and D. 
Dimitrov for discussions. 
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